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Abstract This paper mainly studies nonnegativity decision of forms based on variable substitutions. 
Unlike existing research, the paper regards simplex subdivisions as new perspectives to study variable 
substitutions, gives some subdivisions of the simplex T„, introduces the concept of convergence of the 
subdivision sequence, and presents a sufficient and necessary condition for the convergent self-similar sub- 

Q division sequence. Then the relationships between subdivisions and their corresponding substitutions are 
established. Moreover, it is proven that if the form F is indefinite on T n and the sequence of the successive 
L-substitution sets is convergent, then the sequence of sets {SLS' m ' (F)}^ =1 is negatively terminating, 
and an algorithm for deciding indefinite forms with a counter-example is obtained. Thus, various effective 
substitutions for deciding positive semi-definite forms and indefinite forms are gained, which are beyond 
• the weighted difference substitutions characterized by "difference" . 
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£2 ■ 1 Introduction 

"xT ■ 

Theories and methods of nonnegative polynomials have been widely used in robust control, 
non-linear control and non-convex optimization [H, 0, etc. Some famous research works on 
nonnegativity decision of polynomials without cell-decomposition were given by Polya's Theorem 
q \ 0, [j| and papers [s jz K 

Recently, Yang E| 1> E3| introduced a heuristic method for nonnegativity decision of polynomi- 
.£h ! alSj which is now called Successive Difference Substitution (SDS). It has been applied to prove a 

great many polynomial inequalities with more variables and higher degrees. Yao [Tl| investigated 
the weighted difference substitutions instead of the original difference substitutions, and proved 
that, for a form (namely, a homogeneous polynomial) which is positive definite on K™ , the corre- 
sponding sequence of SDS sets is positively terminating, where R" = {(x±,X2, ...,x n )\xi > 0,i = 
1, 2, ...,n}. That is, we can decide the nonnegativity of a positive definite form by successively 
running SDS finite times. The research results above are all confined to the variable substitutions 
characterized by "difference" . 

Unlike existing research, this paper regards simplex subdivisions as new perspectives to study 
variable substitutions, and obtains various effective variable substitutions for deciding positive 
semi-definite forms and indefinite forms on R™, which are beyond the weighted difference sub- 
stitutions characterized by "difference" . The paper is organized as follows. Section 2 gives some 
subdivisions of the simplex T n , and establishes the relationships between them and their corre- 
sponding substitutions. Section 3 introduces the concept of the termination of {SLS^ m ^(F)}^ =1 , 
which is directly related to the positive semi-definite property of the form F. Section 4 proves 
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that if the form F is indefinite on R™ and the sequence of the successive L-substitution sets is 
convergent, then the sequence of sets {SLS^ m - ) (i 7, )}^ ) =1 is negatively terminating. An algorithm 
for deciding an indefinite form with a counter-example is presented in Section 5, and by using the 
obtained algorithm, several examples are listed in Section 6. 

2 Simplex subdivisions and the corresponding substitutions 

We first introduce some definitions and notations. 

n 

Definition 2.1. Let V = [vij] be an n x n matrix. If ^ Vij = 1, j = 1, 2, • • -n, V is called a 

i=i 

normalized matrix. And the corresponding linear transformation 

is called a normalized substitution, where X,Tg R™, X Tt , T Tr are respectively the transposes 
of X, T, and V is called the substitution matrix of ([T|). 

Lemma 2.1. Let U = [v,ij] = V\V 2 ■ ■ ■ Vfe. If V\, V 2l ■ ■ ■ ,Vk are all normalized matrices, then U is 
a normalized matrix , that is, 

n 

^Uij = l(j = 1,2, • ■■n). 

n n 

And let (xi,x 2 , ■ • • ,x n ) Tl = U(ti,t 2 , ■ ■ ■ , i n ) Tr , then ^ Xj = 1 iff U = 1 - 

i=l i=l 

The proof of Lemma [27X1 is very straightforward and is omitted. 

Definition 2.2. If the form F(X) > for all X e K™ , then F is called positive semi-definite 
on R™ , and the set of all positive semi-definite forms is denoted by PSD; If F(X) > for all 
X(^ 0) G R™ , then F is called positive definite on R™ (briefly, a PD); If there are X and F e R™ , 
such that F(X) > and F(Y) < 0, then F is called indefinite on R^. 

The (n — l)-dimensional simplex is defined as follows 

n 

T„ = {{xi,x 2 , ■ ■ ■ ,x n )\^Xi = 1, (x 1 ,x 2 , • • • , x n ) G R+}. 
t=i 

Definition 2.3. If the form F(X) > for all X G T„, then F is called positive semi-definite 
on T„; If F(X) > for all X G T n , then F is called positive definite on T„; If there are X and 
y G T„ such that F(X) > and F(Y) < 0, then F is called indefinite on T„. 

It is easy to get the following conclusion. 

Lemma 2.2. The form F is positive semi-definite (positive definite, indefinite) on R™ iff F is 
positive semi-definite (positive definite, indefinite) on T„. 

According to Lemma 12.21 for brevity, we suppose the form F is defined on T n in the following 
sections. 

2.1 The barycentric subdivision and the weighted difference substitu- 
tions 

In this subsection, we focus on the barycentric subdivision and its corresponding substitutions. 



Figure 1: the barycentric subdivision 



See Fig.l. The first barycentric subdivision of the simplex T3 consists of 6 subsimplexes. 
Consider the subsimplex T x AO labeled f first, where Ti = (1, 0, 0), A = (|, \, 0), O = (|, |, |) in 
the TiT^T^-coordinate system. Consruct the following substitution 
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By ©, if X = Ti = (1,0,0), then T = (1,0,0); If X = A = (|, |, 0), then T = (0,1,0); If 
X = O = (|, -|, -|) , then T = (0,0, 1). And it is indicated that the subsimplex T\AO and the 
substitution @ correspond to each other. 

Analogously, the other five subsimplexes labeled as 2-6 correspond to the five substitutions, 
respectively, whose substitution matrices are as follows 
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The six substitutions above are just the weighted difference substitutions for n — 3, and the set 
which consists of all the substitutions is called the weighted difference substitution set. Therefore, 
the fisrt barycentric subdivision of T3 corresponds to the weighted difference substitution set, 
that is, there is a one-to-one correspondence between the subsimplexes of the first barycentric 
subdivision of T3 and the weighted difference substitutions. 

Definition 2.4. Let PW3 be the weighted difference substitution matrix set for n = 3. The set 
of linear transformations 



Tr 



B [ai} B [a 2 ] 



B[a m ]T \B[ ai ] ePW 3 }, 



is called the m-times successive weighted difference substitution set, which consists of 6 m substi- 
tutions. 

By Lemma [2. II and Definition ^. 41 we have that the m-th barycentric subdivision of the simplex 
T3 corresponds to the m-times successive weighted difference substitution set. 



2.2 Some more subdivisions and the corresponding substitutions 

Next, we'll present some more subdivisions of the simplex T3 and their corresponding substi- 
tutions. 




Figure 2: Some subdivisions 



Consider Fig.2(a) first. Let T x = (1,0,0), T 2 = (0,1,0), T 3 = (0,0,1), O = (a ,b ,c ), 
A = (ai,6i,0), B = (0,62,02), C = (03,0,03). And the first subdivision of the simplex T 3 
consists of six subsimplexes labeled as 1-6, which correspond to the following substitution matrices, 
respectively. 
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then 



the subdivision is just the barycentric subdivision. 

In Fig. 2(b), the first subdivision of the simplex T 3 consists of four subsimplexes labeled as 1-4, 
which correspond to the following substitution matrices, respectively. 
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In Fig. 2(c), the first subdivision of the simplex T 3 consists of nine subsimplexes labeled as 1-9, 
which correspond to the following substitution matrices, respectively. 
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Next, we'll give the definitions of the self-similar subdivision sequence and the successive sub- 
stitution set. 

Definition 2.5. Given a simplex K. Let sd°{K) = K, and sd l (K) = sd^d^ 1 ^)) , where sd l (K) 
is the subdivision of sd J_1 (i^) for i = 1,2, For an arbitrary subsimplex a of sd l (K)(i = 
0, 1, ••■), if the vertex coordinates of all the subsimplexes of sd(cr) in the a-coordinate system 
equal to the corresponding ones of all the subsimplexes of sd(K ) in the if-coordinate system, then 
sd J (-ftT) is called the z-th self-similar subdivision of the simplex K, and the subdivision sequence 
{sd l (T„)}^ 1 is called self-similar. 



Fig. 3 shows the first and second self-similar subdivisions of the simplex T 3 . 



(a) (b) (c) 



Figure 3: the self-similar subdivisions of the simplex T3 



Definition 2.6. Let PA be the L-substitution matrix set. The set of linear transformations 

{X Tl = A {1] A [2] ■ ■ ■ A [m] T^\Au G PA, i = 1,2, • • • , m}, 
is called the m-times successive L-substitution set, briefly, the m-times SLS set. 

Suppose that the the first subdivision of the simplex T3 corresponds to the normalized L- 
substitution set. By Lamma 12.11 Definition 12.51 and Definition 12.61 we conclude that the m-th 
self-similar subdivision of the simplex T3 corresponds to the m-times successive L-substitution 
set. 

Easily, all the subdivisions of the simplex T3 in the paper can be expanded to the case of the 
(n — l)-dimensional simplex T„. 



2.3 Convergence of the subdivision sequence of a simplex 

In this subsection, we'll introduces the concept of the convergence of the subdivision sequence, 
and presents a sufficient and necessary condition for the convergent self-similar subdivision se- 
quence. 

Definition 2.7. Let a be a subsimplex of T„, the maximum distance between vertexs of a is 
called the diameter of a. 

Definition 2.8. Let K and sd l (K) be defined as Definition 12.51 If for all s > 0, there exists 
N £ N, such that all the diameters of the subsimplexes of sd N (K) are less than e, the subdivision 
sequence {sd^Lf)}?^ is called convergent. 

Definition 2.9. Suppose that the subdivision scheme through which sd l_1 (_ftT) is subdivided 
into sd l (K) corresponds to the substitution set Li for i = 1, 2, • • • . If the subdivision sequence 
{sd l (K)}°^ l is convergent, then the sequence of substitution sets {Li]°^ l is called convergent, and 
if Lj = L, i = 1, 2, • • • , briefly, we say that the sequence of the successive L-substitution sets is 
convergent. 

Next, we'll consider the convergence of the sequence of the successive weighted difference 
substitution sets. 

Lemma 2.3. 1^, HI Let if be a complex. If Kn is the fc-th barycentric subdivision of K, then 
for all e > 0, there exists N € N, such that all the diameters of the subsimplexes of Km are less 
than s. 

By Lemma 12.31 , we have the following theorem. 

Theorem 2.1. The barycentric subdivision sequence {Ki]°^L l of T n and the corresponding se- 
quence of the successive weighted difference substitution sets are convergent. 



Now, we present a sufficient and necessary condition for the convergent self-similar subdivision 
sequence, which plays important roles for nonnegativity decision of forms. 



Theorem 2.2. Let {sd J (T n )}g,j be the self-similar subdivision sequence of T„, then it is conver- 
gent iff an arbitrary f -dimensional proper face of subsimplcxcs of sd(T„) is not the 1-dimensional 
one of T„ . 



Proof. The proof of necessity is straightforward by Definition 12.51 and Definition 12.81 Next, we 
prove sufficiency. Let d be the diameter of T n , and r be the maximum diameter of subsimplexes 
of sd(T n ), then r < d. Since an arbitrary 1-dimensional proper face of subsimplexes of sd(T n ) is 
not the 1-dimensional one of T„, we have r < d. Note I = then I < 1. Obviously, the maximum 
diameter of subsimplexes of the iV-th self-similar subdivision of T n is l N d. When N — > oo, then 
I d — ► 0, so the self-similar subdivision sequence {sd 8 (T n )}^. 1 is convergent. □ 

Consider the following subdivisions of T3 in Fig. 4. 

In Fig. 4(a), the first subdivision of the simplex T3 consists of 3 subsimplexes labeled as 1-3, 
which correspond to the following substitution matrices, respectively. 



(8) 



In Fig. 4(b), the first subdivision of the simplex T3 consists of 5 subsimplexes labeled as 1-5, 



which correspond to the following substitution matrices, respectively. 
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In Fig. 4, it is straightforward that the 1-dimensional proper face T\Ti of the subsimplex T1T2O 
is just the one of T1T2T3. By Theorem l2.21 we have that the two self-similar subdivision sequences 
of T n aren't convergent. 





(a) (b) 

Figure 4: Some more subdivisions 



3 Termination of the Sequence of the Successive Substitu- 
tion Sets of a Form 

In this section , we'll define the termination of the sequence of the successive substitution sets 
of a form, which is directly related to the positive semi-definite property of the form F. 
Let PA = {Arq \i = 1, 2, • • • , k} be the L-substitution matrix set. 

Definition 3.1. Given the form F on T„, when [ai], [a 2 ], • • • , [a m ] traverse all the memebers of 
1, 2, • • • , k respectively, we define the set 

SLS^(F)= |J ■■■[J[J% ll yl hl ...A k| I T ). 

lat m ] [aj] [01] 

which is called the m-times successive L-substitution set of the form F. 



Definition 3.2. We define the sequence of sets {SLS (m) (F)}^ =1 as follows 

{SLS(F) (m) }~ =1 = SLS(F),SLS (2) (F),--- . 

Let a = (oti, a 2 , • ■ • , a n ) e N n , and let |a| = ai + a.% + • • ■ a n . Then we write a form F with 
degree d as 

r — 2_ t c a x i %2 ' ' ' X n ■ 
\a\=d 

Definition 3.3. The form F is called trivially positive if the coefficients c a of every term 
x^x^ 2 ■ ■ ■ x"" in F are nonnegative; If -,•••,-)< 0, then F is called trivially negative. 

Lemma 3.1. Given the form F on T n , if the form F is trivially positive, then F £ PSD; If the 
form F is trivially negative, then F € PSD. 

Definition 3.4. Given a form F on T„, if there is a positive integer k such that every element 
of the set SLS (fe) (F) is trivially positive, then the sequence of sets {SLS (m) (-F)}m=i is called 
positively terminating; If there is a positive integer k and a form G such that G € SLS^ (F) and 
G is trivially negative, the sequence of sets {SLS'™' (F)}^ =1 is called negatively terminating; The 
sequence of sets {SLS*-" 1 ' {F)}^ =1 is neither positively terminating nor negatively terminating, 
then it is called not terminating. 



4 Nonnegativity decision of forms 

Given a form F on T n . We know that if the sequence of sets {SLS*- m '(i 7 ')}^ =1 is negatively 
terminating, then we can conclude that the form F is positive semi-definite. Thus there is a natural 
question, that is, which kind of forms can be solved by the method? The following theorem answers 
the question. 

Theorem 4.1. Let the form F be positive definite on T„. If the sequence of the successive 
^-substitution sets is convergent, then the sequence of sets {SLS^™' (i r )}^_ 1 is negatively termi- 
nating. 

Proof. We only give the proof for the ternary form with degree d, and the multivariate form can 
be gotten by induction. 

Suppose that F{x\, x 2 , x 3 ) — aijkx\x 3 2 x 3 . An arbitrary m-times successive L-substitution 

i+j+k—d 

can be written as 

X\ = fcuui + ki 2 u 2 + k 13 u 3 , 

x 2 = k 2 \ui + k 2 2U 2 + k 23 u 3 , (10) 
x 3 = k 31 ui + k 32 u 2 + k 33 u 3 . 

Let 

#1 = &U +&21 + &31, K 2 = k 12 + k 2 2 + k 32 , K 3 = k 13 + k 23 + k 33 , 

h = Kiui, t 2 =K 2 u 2 , t 3 = K 3 u 3 , 

Kl — Kl , 2 if ! ' 3 fc 1 , 

Qi = |f-fei, a 2 = ^-k 2 , a 3 = ^-k 3 , 

Pi = ^-h, P 2 = ^-k 2 , Ps = ^-k 3 . 

then (JTUJ) becomes 



kih + (ki + ai)t 2 + (ki + [3i)t 3 , 

k 2 h + (k 2 + a 2 )t 2 + (k 2 + /3 2 )t 3 , (11) 
Mi + (k 3 + a 3 )t 2 + (k 3 + (3 3 )t 3 , 



3 3 3 

where fc, = 1, ^ a, = and ^ — 0. 

i— 1 i— 1 z— 1 



Let t = ti + 12 + £3, pip can be rewritten as 

x\ = kit + a 1 t 2 + 0it 3 , 

x 2 = k 2 t + a 2 t 2 + 2 t3, (12) 
£3 = M + a 3 t 2 + 3 t 3 . 

Thus 

h, t 3 ) = F{k x t + a x t 2 + 0it 3 , k 2 t + a 2 t 2 + 2 t 3 , k 3 t + a 3 t 2 + 3 t 3 ) 

= E ai k{kit + ai t 2 + 0it 3 )\k 2 t + a 2 t 2 + hh) 1 {k 3 t + a 3 t 2 + 3 t 3 ) k 

= E a ^(k\f+ E zjij*?"!^*!^)- 

•i+j+fc— d p+g+r— i,p^i 

+ E ^^!/WD(^' + E ^ fc 3^ft r * ? ^D 

= E a ijk k\k 3 2 k%t d + E (/> l3 k(ki,k2,k 3 ,a 1 ,a 2 ,a 3 ,f3 1 ,l32,f3 3 )t\t : 2 t 3 : 

i-\-j-\-k—d i J rj-\-k=d 

= F(ki, k 2 , k 3 )t d + E <l>ijk(fa,fo,k3,ai,oi2,a3,l3i,fc,fh)t\1%t$ 
= E ( ^nTl^^ 1 '^ 2 '^ 3 ^ + rij vk{ki 1 k 2 ,k 3l ai,a 2 ,a 3l Pi 1 2 ,0 3 )}t\t : 2 t 3 



A\j\k\ 

i-\-j-\-k—d 
i-\-j-\-k—d 



where 



(13) 



= VmTj F(fci, fe) + (j>ijk{kl, k 2 , k 3 , Oil, «2, 03,01,02, 03). (14) 
Z.J 

Obviously, 

lim <pijk(ki,k2,k3,ai,a2,a 3 ,0i,02,03) = 0. (15) 

And since -F(xi, X2, X3) is positive definite on T„, there exists e > 0, such that 

F(fe 1 ,ft 3) * 3 )>e>0. (16) 

On the one hand, the vertexs of the subsimplex which corresponds to the successive L- 
substitution (fTTj) are respectively 

(ki,k 2 , k 3 ), (ki + ai,k 2 + a 2l k 3 + a 3 ), {k x + 0i,k 2 + 02, k 3 + 3 ). 

By Theorem (|2.1|) . OLi, a2,ot3,0i,0 2 ,03 can be sufficiently small when m is sufficiently large. 

On the other hand, for F is continuous on T n and by (jT4j) - (jTB]) , we have > when 

tti, 0(2, as, 0i, 02, 03 are sufficiently small. 

Putting together the above two aspects, we have that there exists a sufficiently large integer m, 
such that F becomes trivially positive by (|11[) . For the successive L-substitution (flU)) is arbitrary, 
the sequence of sets {SLS*-" 1 ' (F)}^ =1 is positively terminating. □ 

Le t L be the Yang-Yao's substitution set in Theorem 14. 11 then the theorem is the main result 

in 



ll|. However, the proof of the theorem in this paper is different from the one given in 11 1. 



According to the proof of Theorem 14. 11 we obtain the following conclusion. 



Corollary 4.1. Let the form F be positive definite on T„. If the sequence of the successive 
L-substitution sets is convergent, then by an arbitrary m-times successive L-substitution, when 
m is sufficiently larger, F can become a nonlacunary form whose coefficients are all positive. 

In Theorem 14.11 if the sequence of the successive L-substitution sets isn't convergent, the 
conclusion of the theorem is not always true. 

Theorem 4.2. Let the form F be positive definite on T n . If the sequence of the successive L- 
substitution sets isn't convergent, the sequence of sets {SLS (m) (F)}^ =1 isn't always negatively 
terminating. 

Proof. Let L be the substitution set, which correspands to the subdivision of the simplex T3 given 
by Fig. 4(a). And we have concluded that the sequence of the successive L-substitution sets isn't 
convergent. 

Given the form 

F(x 1 ,x 2 ,x 3 ) = (xi - x 2 + x 3 ) 2 + x%, (x 1 ,x 2 ,x 3 ) e T 3 . (17) 
Apparently, F is positive definite on T„. Consider the m-times successive L-substitution 
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And without difficulty, we have 
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Then 



F(Afa ■ {x 1: x 2 ,x 3 ) T ) =x\ - 2 Xl x 2 + 7^x lX3 + 2x 2 2 + (1 

1 5 2_ 2 
4^ 32m 3m + 1 ) x 3- 



—)x 2 x 3 + 



(20) 



When m — > 00, the coefficient of the term x\x 2 is -2. Therefore, the sequence of sets {SLS*-™-* (F)} 



isn't terminating. 



□ 



5 Decision of indefinite forms 

Many problems, such as the inequality disproving, are always transformed into decision of 
indefinite forms. Given a form F on T n . Suppose that there exists Xq € T n such that F(Xq) > 0. 
It is well-known to us that if the sequence of sets {SLS^^L)}^]^ is negatively terminating, then 
the form F is indefinite. Then it follows a question naturally: for a indefinite form f on T„, is 
the {SLS*- m ' ) (i 7, )}^ =1 negatively terminating? The question is answered by the following theorem. 

Theorem 5.1. Let the form F be indefinite on T„. If the sequence of the successive L-substitution 
sets is convergent, then the sequence of sets {SLS^ m ^(i r )}^_ 1 is negatively terminating. 



Proof. Let PA be the substitution matrix set which corresponds to the substitution set L. Since 
the form F is indefinite on T„, there exists Xq G T n such that F(Xo) < . And F is continuous 
on T n , so there exists a neighborhood U(Xq) C T n of Xq (If Xq is on the boundary of T n , then 
we take U(X ) H^") sucn tnat < for all X G U(A ). For the sequence of the successive 

L-substitution sets is convergent, then there exists a subsimplex er of the fc-th subdivision of T„, 
which corresponds to the fc-times successive L-substitution 

XTl = B [ti] B [i2] ■ ■ ■ B [ik] TTl , B [ii]i B [i^ ■ • • i B [i k ] G PH7 «' 

satisfying a C U(Ao), where A; is a sufficiently larger integer. Thus, —F(X) is positive definite 
on a. By Theorem 14. 11 the sequence of sets {SLS'-" 1 -' (— F)}^ =1 is positively terminating, so there 
exists a Ltimes successive L-substitution 

* Tl ' = B M B M ■ ■ ■ B {jl] T^, B W , B [J2] , • • ■ , B m G PW n , 

satisfying that i ?1 (.Brj 1 ].Bu 2 i • ■ ■ Br^iSy lSy-i • • ■ SujT ) is trivially negative. Therefore, the se- 
quence of sets {SLS^ m '(F)}^ =1 is negatively terminating. □ 

By the proving process of Theorem (|5.ip . we obtain the following algorithm, which is used to 
decide the indefinite form with a counter-example. 
Algorithm (SLS) 

Input: the form F G Q[xi,X2, • • • , x n ], where F is positive definite or indefinite on R" . 

Output: "F GPSD", or U F(X ) < 0". 

stepl: Let F = {F}. 

step2: Compute (J SLS(F), 

FeF 

Let 

F = (J SLS(F) - {trivially positive forms in (J SLS(F)} = {F {1] , F [2] , •••,%}, 

FeF FeF 

where 

F M = F(S H X Tr ), S M G PW„. 

step3: Let 1 = [[1],[2],--- ,[k]]. 



step31: If F is null, then output "F ePSD", and terminate. 
step32: If there is a trivially negative form F^ G F, then output u F(Xo) < 0", 
and terminate, where 

= B m(~' -> -) Tr > B m = B mm] B imm] • ■•- B [-fMH]> 

is the z-th component of /, is the j-th component of and m is the 

the total number of the components of 

step33: Else, Compute [j SLS(F). Let 

FeF 

F = (J SLS(-F) - {trivially positive forms in |J SLS(F)} 

FeF FeF 

= {P[op(/[l]),l],- ■ ' ,i ? [op(/[l]),i 1 ],-P , [op(/[2]),l], ' ■ ' ,P[op(/[2]),i 2 ], 
' ' ■ ,P[op(J[fc]),l],' ' ■ ,P[op(J[fc]),i fc ]}, 

where 

•^[opC/W)^] = F ( B [o P (L[i\)j]X Tr ), 
and op(J[i]) extracts operands from op(J[i]). And let 



L = [[op(/[l]), 1], ■ ■ • , [op(/[l]), ij, [o P (/[2]), 1], • • ■ , [op(/[2]),i a ], 
••• ,[op(I[fc]), !],••• ,[op(I[k])M, 



then go to step3. 

By Algorithm SLS, we design a Maple program called SLS, see Appendix. To the program 
SLS, there are some positive semi-definite forms making the program do not terminate, that is, 
we cann't decide these positive semi-definite forms by the method. 

6 Examples 

In this section, we demonstrate the program SLS with some examples on a computer with 
Intel(R) Core(TM)2 Duo CPU (E7200 @ 2.53GHz) and 3.25G RAM. 

Let Li,L 2 ,L 3 be the substitution sets which correspond to the subdivisions of the simplex T 3 
given by Fig.l, Fig. 2(b) and Fig. 2(c), respectively. And let PAi be the Li-substitution matrix set 
for i = 1,2,3. 

Example 1. Show that the following form is positive scmi-dcfinitc on 

F(x, y, z) = x(x - y) 5 - y(-z - y) 5 - z(x - zf. (21) 

Utilize the program SLS and execute order SLS(.F, PAi, [x, y, z]) for i = 1,2, 3. And the running 
results are shown in Table 1. So F positive semi-definite on 



The substitution set 


Running SLS times 


Output 


CPU time (s) 


Li 


3 


F ePSD 


0.016 


L 2 


3 


F ePSD 


0.015 


L 3 


2 


F ePSD 


0.032 



Table 1. 



Example 2. Show that the following form is indefinite on 

F(x, y, z) = x 4 y 2 - 2x 4 yz + x 4 z 2 + 3x 3 y 2 z - 2x 3 yz 2 - 2x 2 y 4 - 2x 2 y 3 z + x 2 y 2 z 2 + y 6 , (22) 

Execute order SLS(.F, PAi, [x, y, z]) for i = 1, 2, 3, and the running results are shown in Table 
2. Obviously, F(0, 1, 0) = 1 > 0, so the form F is indefinite on R\. 



The substitution set 


Running SLS times 


Output 


CPU time (s) 


Li 


3 


F (37 91 85 \ - n 
r ^81 ' 324' 324 J ^ u 


0.094 


L 2 


3 


-^(24' 3' 24) ^ ^ 


0.062 


L 3 


2 


^(27' 27' 27) ^ ^ 


0.094 



Table 2. 



Example 3. Prove or disprove that, for x > 0, y > 0, z > 0, and x + y + z ^ 0, 

2, x 2 y 2 z 2 s ,x e +y 6 + z 6 ,i n 

77 + — — + - ~, F >0, (23) 

3 y + z z + x x + y' v 3 ; v ; 

Take off denominators of the left polynomial, and denote the new polynomial by F. Execute 
order SLS(F, PAi, [x, y, z\) for i = 1,2,3, the running results are shown in Table 3. So the 
inequality doesn't hold. 



The substitution set 


Running SLS times 


Output 


CPU time (s) 


Li 


5 


EV2159 3685 3661 \ „ n 
r V5832' 11664' 11664 ^ ^ u 


582.750 


L 2 


5 


t?( 7 37 31 \ ^ n 
^ V 24 ' 96 ' 96 > ^ u 


21.469 


L 3 


3 


P(31 25 25 \ ^ n 

r y&i ' si ' si 1 ^ u 


234.484 



Table 3. 



The examples above indicate the effectiveness of Algorithm SLS. Thus, we obtain various 
effective substitutions for deciding positive semi-definite forms and indefinite forms which are 
beyond Yang's substitutions characterized by "difference" . 
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Appendix. Maple Program SLS 



SLS : =proc (poly , A , var) 

local f,i,j,k,m,n,p,s,t ,F,G,H,M,newvar , st , Var : 
uses combinat, Linear Algebra : 
t:=time() : F:=[[poly, [0]]] :n:=nops(var) :r=100 
Var :=convert( var, Vector) : 
for i to nops(A) do 
for j to n do 

newvar [i , j] :=op( j , convert (A [i] .Var, list) ) : 

od: 

od: 

for s to r do 

m:=nops(F): f: = []: 
for k to m do 
G[k]: = []: 

for i to nops(A) do 

st:={seq(Var[j]=newvar[i, j] ,j=l. .n)>: 

G[k] : = [op(G[k]), [expand(subs(st,F[k] [1] ) ) , [op(F [k] [2]),i]]] : 

od: 

od: 

F: = [seq(op(G[u]) ,u=l. .nops(F))] : 
for i to nops(F) do 

if max( [coeff s(F[i] [1] )] ) < then 
print (F[i] [2]) : 
M :=IdentityMatrix(n) : 
for j from 2 to nops (F [i] [2] ) do 
M:=M.A[F[i] [2] [j]] : 

od: 

print (convert (M. Vector [column] (n, 1/n) ,list)) : 
print (time ()-t , ' second' ) : 
return ("The form is indefinite") : 
elif min( [coeffs(F[i] [1] )] ) >=0 then 
f : = [op(f),i] : 

fi: 

od: 

if nops(f)>0 then 

F:=subs({seq(F[f [p] ]=NULL,p=l . .nops (f ) ) } ,F) : 

fi: 

if nops(F)=0 then 

print(s) : print (time ()-t, 'second') : 
returnO'The form is positive semi-definite") : 

fi: 

od: 

end proc : 



